{\rtf1\ansi\ansicpg1252\cocoartf1344\cocoasubrtf720
{\fonttbl\f0\fswiss\fcharset0 Helvetica;}
{\colortbl;\red255\green255\blue255;}
\paperw11900\paperh16840\margl1440\margr1440\vieww15960\viewh14820\viewkind0
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\f0\fs36 \cf0 Exerc\'edcios da Aula 4 - An\'e1lise explorat\'f3ria de dados
\fs24 \
\

\fs36 Exerc\'edcio 4.2\

\fs24 \
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs28 \cf0 cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")\
cervejas=as.factor(cervejas)\
cervejas1=(sort(table(cervejas), decreasing=F))\
#class(cervejas1)\
par(mfrow=c(1,2))\
# Representa\'e7\'e3o - Gr\'e1fico de barras\
barplot(cervejas2, xlab= "Cervejas", ylab="Frequ\'eancia")\
# Representa\'e7\'e3o - "Dotplot"\
dotchart(cervejas2, xlab="Frequ\'eancias")\
par(mfrow=c(1,1))\
# O gr\'e1fico "dotplot" apresenta maior raz\'e3o dado/tinta.
\fs24 \
\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs36 \cf0 Exerc\'edcio 4.3\
\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs28 \cf0 caixeta <- read.csv("caixeta.csv", as.is=T)\
caixeta$dap=caixeta$cap/pi\
names(caixeta)\
# Criando o histograma dos fustes dos caixetais\
hist(caixeta$dap, main = "Histograma de DAP", xlab = "Dap", ylab = "Fuste")\
# Criando o histograma da altura para os diferentes caixetais\
library(lattice)\
histogram( fuste ~ dap | local , data=caixeta, main = "DAP por local", xlab = "Dap", ylab= "Fuste")\
\
# Sim, h\'e1 diferen\'e7as nas estruturas dos DAPs dos fustes entre os tr\'eas locais, visualmente \'e9 \
# poss\'edvel inferir que h\'e1 mais fustes com menor DAP no Retiro e em Chauas em rela\'e7\'e3o \'e0 Jureia.\
\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs36 \cf0 Exerc\'edcio 4.4\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs28 \cf0 \
egrandis <- read.csv("egrandis.csv", sep = ";", as.is=T)\
class\
#head(egrandis)\
#names(egrandis)\
# boxplot dos DAPs por regi\'e3o\
boxplot(egrandis$dap~egrandis$regiao)\
# boxplot dos DAPs por rota\'e7\'e3o\
boxplot(egrandis$dap~egrandis$rotacao)\
\
#mean(egrandis$ht)\
# Avaliando a normalidade\
qqnorm(egrandis$ht, main = "Norm alturas (quantil-quantil)", xlab = "Quantis te\'f3ricos", ylab = "Quantis amostrados")\
qqline(egrandis$ht)\
\
# Ao analisar o gr\'e1fico quantil-quantil, visualmente, a altura n\'e3o segue uma distribui\'e7\'e3o normal.\
\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs36 \cf0 Exerc\'edcio 4.5\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs28 \cf0 caixeta <- read.csv("caixeta.csv", as.is=T)\
caixeta$dap=caixeta$cap*pi\
caixeta.Tab=caixeta[caixeta$especie=="Tabebuia cassinoides", ] # Selecionar somente as Tabebuia\
\
## Agrupar os fustes por \'e1rvore\
caixeta.Tab$ab=pi*((caixeta.Tab$dap/2)^2)\
caixeta.arv.Tab= aggregate(caixeta.Tab$ab,by=list(caixeta.Tab$arvore,caixeta.Tab$local,caixeta.Tab$h),FUN=sum)\
colnames (caixeta.arv.Tab) = c("arvore", "local", "h", "ab")\
#head(caixeta.arv.Tab)\
#tail(caixeta.arv.Tab)\
caixeta.arv.Tab$dap=2*sqrt(caixeta.arv$ab/pi)\
\
## Plotar a rela\'e7\'e3o altura x \'e1rea basal x local\
coplot( h ~ dap | local , data=caixeta.arv.Tab, main = "Rela\'e7\'e3o DAP x Alturua / Local", xlab = "Altura", ylab = "DAP")\
#Pelo library\
library(lattice)\
xyplot( h ~ dap | local , data=caixeta[caixeta$especie=="Tabebuia cassinoides", ], main = "Rela\'e7\'e3o DAP x Alturua / Local", xlab = "Altura", ylab = "DAP")\
\
## Verificando a linearidade - Indexando dados\
caixeta.arv.Tab.chauas=caixeta.arv.Tab[caixeta.arv.Tab$local=="chauas", ]\
caixeta.arv.Tab.retiro=caixeta.arv.Tab[caixeta.arv.Tab$local=="retiro", ]\
caixeta.arv.Tab.jureia=caixeta.arv.Tab[caixeta.arv.Tab$local=="retiro", ]\
# Verificando a linearidade - Gr\'e1ficos\
par(mfrow=c(1,3))\
scatter.smooth(caixeta.arv.Tab.chauas$dap, caixeta.arv.Tab.chauas$h, main = "Linearidade Chaua", xlab = "DAP", ylab="Altura")\
scatter.smooth(caixeta.arv.Tab.retiro$dap, caixeta.arv.Tab.retiro$h, main = "Linearidade Retiro", xlab = "DAP", ylab="Altura")\
scatter.smooth(caixeta.arv.Tab.jureia$dap, caixeta.arv.Tab.jureia$h, main = "Linearidade Jur\'e9ia", xlab = "DAP", ylab="Altura")\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs36 \cf0 \
\
\
\pard\tx566\tx1133\tx1700\tx2267\tx2834\tx3401\tx3968\tx4535\tx5102\tx5669\tx6236\tx6803\pardirnatural

\fs28 \cf0 \
}